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Abstract 

A mixture of the multicomponent Bose-Einstein condensate is considered, where 
each component moves with its own velocity. As a result of the relative motion, the 
mixture stratifies when the relative velocity reaches a critical value. Stability conditions 
for a binary moving mixture are derived and the critical velocity is found. 
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Ultracold atomic gases, forming Bose-Einstein condensates, arc nowadays an active field 
of research (see reviews [1-5]). Dilute Bose gas can be trapped and cooled down to temper- 
atures when practically all atoms are in Bose-Einstein condensate. One can create mixtures 
of several Bose-condensed species. The stability of such mixtures depends on the mutual 
atomic interactions. The condition of stability for binary mixtures was, first, found by 
Nepomnyashchy [6], considering the collective excitation spectrum in the Bogolubov and 
long-wave approximations. When this stability condition is not valid, the mixture is not sta- 
ble and stratifies onto separate components. The analysis of Nepomnyashchy [6] was based 
on the Green functions approach for a superfiuid mixture with broken gauge symmetry. As 
was shown later [7] , the same stability condition holds true for a fluid mixture with conserved 
gauge symmetry. 

Recently, it has become feasible to create moving wave packets of Bose-Einstein conden- 
sates, which could collide and move through each other [8-12]. Therefore, it is interesting to 
study the interplay between the possibility of different Bose components to mix and the in- 
fluence on their stability of their relative motion through each other. Such a relative motion 
of two fluids could change the stability condition of mixing [6,7]. Bose-Einstein condensates, 
formed in modern traps, are made of dilute gases, for which the Gross-Pitaevskii equation 
provides a very good description [1-5]. Actually, this equation is an exact equation for the 
coherent system at zero temperature [13]. The aim of the present communication is to an- 
alyze a mixture of Bose-Einstein condensates, whose components move through each other, 
under conditions of low temperature, when the system is described by the Gross-Pitaevskii 
equation. This situation corresponds to current experiments with Bose-Einstein condensates 
of dilute gases. 

An atomic gas, of N atoms in volume V, is dilute when its density p = N/V and the 
s-wave scattering length a are connected by the inequality pa^ <^ 1. This can be generalized 
to a mixture of species enumerated by the index i — 1,2, ... as the inequality ^piPjO^^ <C 1, 
in which an = a.i and are the corresponding scattering lengths, and Pi = Ni/V. The 
Gross-Pitaevskii equation for a multicomponent mixture is 

d 

Ol Vji^, i) = Hj[ri] rij{r, t) , (1) 
where r]i{v,t) is the coherent fleld of the i-component; in the nonlinear Hamiltonian 

HM = ^ + U, + Y.Mr}3? (2) 

* 3 

the flrst term is the kinetic-energy operator, with p = —ihV, and being atomic mass, C/j 
is an external (trapping) potential, and the third term describes contact atomic interactions 
with 

rriimj 

The normalization condition for each component reads 

\\r)i\\^ = I \r)i{r,t)\'' dr ^ Ni . 
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Elementary collective excitations of the system can be obtained by linearizing the Gross- 
Pitaevskii equation (1) around a given stationary solution. Following here this way, we 
should keep in mind that our aim here is to study a mixture of moving components. Assume 
that each component moves with a fixed velocity Sj. Then taking account of this motion can 
be done by means of the Galilean transformation replacing r)j by rjj exp{^mjSj ■ r}. To find 
collective excitations in the moving mixture, we write the coherent field rjj in the form 

rtj{r, t) = [(^,(r) + u,(r)e— * + v;{r)e^-'] exp (m, s, • r - E,t)^ , (3) 

where (pj{r) satisfies the stationary equation 

+ p • Si + ^ <^i(r) = Ei^i{T) , (4) 



with the normalization Hv^ilP = Ni. Substituting Eq. (3) into the Gross-Pitaevskii equation 
(1) and linearizing it with respect to the functions %(r) and Vj{v), we obtain 

fiuj - p ■ Si - H[ip\ ^ + Eij Ui-ipiY^ Aij [ip*Uj + (fjVjj = , 

(2 \ 
hu}-p-Si + H[ip] + - Eij Vi + ip* Aj + (fjVj) = . (5) 

This is a generalization of the Bogolubov - de-Gennes equations to the case of a mixture of 
moving components. 

Being interested in collective excitations above the ground state of trapped atoms, we 
keep in mind that the wave functions, defined by Eq. (4), are real- valued. Then we can 



resort to the local-density approximation, setting (pi{r) ~ ypj(r) and assuming that Ui and 
Vi are proportional to e*^"". Supposing that atomic interactions are sufficiently strong, so 
that mean potential energies are essentially larger that mean kinetic energies, we may omit 
the terms containing the spatial differentiation of Pi{r). In this approximation. 



+ + (6) 



Then Eqs. (5) reduce to 



?i ( a; - k • Sj - ^ I Ui - YAij./p~p]{uj + vj) = , 



^ l^u; - k • Si Vi^Y Aj^/PiFA^o + "^i) = ■ (7) 

The latter equations possess a nontrivial solution if the related determinant is zero. To write 
down the resulting equation in a compact form, we introduce the notation 

^Nc^^^+f^V> c\^P^Ai^ (8) 



2m,- / m 
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where Ai = An. If there would be no interactions between different species, then each 
component would have the Bogolubov spectrum uji{k) given by Eq. (8). Note that the 
existence of the Bogolubov spectrum for a single-component Bose-Einstein condensate of 
trapped atoms is well confirmed by experiment [14]. Define also 

/ \ 1/2 

a;.,^Q,fc, 4-(;^J (9) 

And let us consider the case of a binary mixture. The spectrum of collective excitations, 
following from Eq. (7), is given by the equation 

^ uj - Elf - ujI\ [(a; - £2)^ - 1^1] = i^u , W 
in which = k ■ Sj. 

It is convenient to connect the system of coordinates with one of the components, say 
with the first one, so that Si = 0. Then S2 = s is the velocity of the second component with 
respect to the first component. In that case, the spectral equation (10) yields 

u"^ - 2eu;^ - (u;f + u;^ - cu^ + 2ziJ\uj + iJ\ (J\ - s^) ^ (11) 

where 

£ = k • s = /cscos'j? . (12) 

Equation (11) has four solutions for oj — a;(k). But only two of them can be stable, being 
positively defined for all ■& e [0, 27r]. The necessary condition for Eq. (11) to possess two 
positive solutions for arbitrary e e [— /cs, ks\ by Descartes theorem, is 

— 03^2 > ujfk'^s'^ . (13) 

In the long-wave limit, when k ^ 0, we have the stability condition 

C1C2 - C12 > . (14) 

This shows that, when the immovable mixture is stable, satisfying the inequality cfcl > 0^2, 
it will retain its stability only up to the critical velocity 



= -^dcl - c\, . (15) 

Cl 



Reaching this relative velocity of one component through another, the binary mixture be- 
comes unstable and stratifies onto two spatially separated components. The stratification 
starts inside the cone of the angle 

■&C — arccos — . (16) 
The stability condition (14), with Eqs. (8) and (9), transforms to 

A1A2 - > — A^s" . (17) 

P2 
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Respectively, the critical velocity (15) can be written as 



Taking into account the relation between Aij and the scattering length a^, the stability 
condition (17) can be presented in the form 

0,10,2 (mi+m2)^ m^ois^ 

— T~ ^ A r A 1-2 9~ • U^J 

ai2 4mim2 47rft ^2^12 

Finally, for colliding wave packets with nonuniform densities Pi(r), the stratification would 
occur not everywhere in space, where the moving packets overlap, but only there, where the 
inequality 

becomes invalid, provided that A1A2 > A^- 

The stability conditions for the mixture of components with relative motion have been 
obtained here by considering the stabihty of the spectrum of collective excitations. This is 
what is called the dynamic stability [2] . There exists as well the notion of thermodynamic 
stability, when one compares the corresponding free energies, or, at zero temperature, the 
related internal energies. For immovable mixtures, these two types of stability result in 
equivalent stability conditions. However, for the mixtures of moving components, these two 
kinds of stability are not equivalent. To show the latter, we may, for simplicity, employ the 
uniform approximation, writing the internal energy of the mixture as 

Emix = ^ (a^i + N^A2 + 2N,N2Au) + K , (21) 

where V — V1 + V2 and the last term is the kinetic energy 

{Nxmxsl + A^2m2S^) . (22) 

If the components are spatially separated, that is, stratified, their total internal energy is 

2Vi 2V2 ^ ' 

Because of the additivity of kinetic energy, the comparison of E'^ix and Egtr does not in- 
volve the velocities of the moving components, so that the thermodynamic stability condi- 
tion Emix < Est,- does not depend on the relative motion of these components. Moreover, 
employing the thermodynamic stability condition, one commonly invokes the condition of 
mechanical equilibrium, according to which the pressures of the stratified components, 

dEstr _ NjA, 
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are equal, so that pi — p2- The latter yields the equality 
from where, recalling that V — Vi + V2, we have the relations 



NU2 Z ^ 1 + 



This makes it possible to present Eq. (23) as 



Estr = ^ {^N^A, + N^A2 + 2N,N2^/A^^ + K . (24) 

Now, comparing Eqs. (21) and (24), one notices that Emix < Egtr if and only if \/ A1A2 > A12. 
However, the usage of the condition of mechanical equilibrium for moving separated com- 
ponents does not look appropriate. The thermodynamic stability condition can be applied 
only to immovable mixtures. 

The method of studying the stability of mixtures with a relative motion of their com- 
ponents by means of the dynamic stability condition, requiring the positive definiteness 
for the spectrum of collective excitations, is quite general and can be applied to arbitrary 
mixtures. In particular, this can also be an atomic-molecular mixture. Such mixtures are 
formed by means of the Feshbach resonance technique [2,15]. Two Bose atoms can form a 
quasi-molecule. Then, it is possible to create a two-component mixture consisting of atoms 
themselves and of the molecules formed by these atoms. When both atoms and molecules 
are in the Bose-condensed state, they again can be described by the system of the Gross- 
Pitaevskii equations for the corresponding coherent fields. Thus, if the interaction amplitude 
for the reaction "atom-|-atom molecule" is denoted by B, then the related Gross-Pitaevskii 
equations are 

d - (9-1 

i^-Q^Vi = Hi['n]rii + Brilr]2 , ih ri2 = H2['n]r]2 + - Brjl , (25) 

where Hi[r]] are the same Hamiltonians (2). When atoms as well as molecules are immovable, 
their mixture is stable under the condition 

B 



A1A2 > A12 + 

where p2 = N2/V . But if one is able to organize the relative motion of the atomic and 
molecular components with respect to each other, then, to find the stability condition for 
such a moving mixture, one has to follow the same procedure as described above. That is, the 
form (3) of the coherent fields has to be substituted into Eqs. (25), after which the spectrum 
of collective excitations can be obtained. Requiring the positiveness of this spectrum yields 
the stability condition under the relative motion of components. 

Concluding, if an immovable mixture of Bose-condensed species is stable, it can stratify if 
the components move with respect to each other. The stratification occurs when the relative 
velocity of motion reaches a critical value. The conditions of stability and, respectively, 
stratification can be derived by linearizing the Gross-Pitaevskii equations for the components, 
with the inclusion of the Galilean transformation for the coherent fields, as is done in Eq. 
(3). 
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